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CONNECTION ADMISSION CONTROL BASED ON 

BANDWIDTH AND BUFFER USAGE 

BACKGROUND 

1. FIELD OF THE INVENTION 

S The present invention pertains to telecommunications, and particular to 

connection admission control aspects of telecommunications traffic management. 

2. RELATED ART AND OTHER CONSIDERATIONS 

In telecommunications, traffic management is the art of providing users with the 
service they need and have paid for. One aspect of traffic management, connection 
10 admission control (CAC), checks to ensure that resource consumption of new 
connections will not violate the quality of service (QOS) requirements of new and 
existing connections before admitting the new connections on the network. Relevant 
resources involved in connection admission control (CAC) are channel numbers, 
bandwidth, and buffer space. 

15 Within connection admission control (CAC), effective/equivalent bandwidth 

methods are based on the asymptotic behavior of a tail of a queue length distribution. 
The algorithms of such methods calculate an effective bandwidth based on the traffic 
descriptor, QOS requirements, and buffer resources. The effective bandwidth methods 
are often restricted to a certain type of traffic sources. Some results apply to on-off 

20 Markov fluids {see Kesidis, Walrand, Chang, "Effective Bandwidths for Multicall 

Markov Fluids and Other ATM Sources", IEEE Trans. Networkings Vol. 1, No. 4, pp. 
424 - 428, August 1993), others to leaky-bucket shaped traffic, see Lo Presti, Zhang, 
"Source Time Scale and Optimal Buffer/Bandwidth Trade-off for Regulated Traffic in 
an ATM Node", UMASS CMPSCI Technical Report UM-CS-96-38, 

25 Some real sources may be closest to the Markov fluid-type model, for instance 

encoded speech. This type of source needs no further shaping to be suitable for transfer 
through a network, although it has no maximum burst size. On the other hand, packet 
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traffic is typically leaky bucket shaped with a peak rate, a sustainable rate, and a 
maximum burst size. The asymptotic behavior of a queue loaded with leaky bucket 
traffic differs significantly from the behavior of a queue loaded with Markov traffic. 
The leaky bucket-loaded queue has an upper bound on its length, while the Markov- 
5 loaded queue has no upper bound. Using the asymptotic Markov-type results on leaky 
bucket traffic will completely loose the leaky bucket burstiness information. Using the 
leaky bucket results with a Markov source does not work because they rely on the 
linoited burst size. 

What is needed, therefore, and an object of the invention, is a connection 
10 admission control (CAC) technique which is effective in multitudinous traffic scenarios. 



BRIEF SUMMARY OF THE INVENTION 

The present invention approximates probabihty of loss using a log moment 
generating function and its two partial derivatives of workload on a queue over a time 
interval. The approximation uses four state variables, which depend on the log moment 
15 generating function and its two partial derivatives. The four state variables are as 
follows: 

z(s,t) Linear term in approximation to log loss ratio at working point (s,t) 

c(s,t) Argument of logarithmic term in approximation to log loss ratio at 
working point (s,t) 

20 B(s,t) Buffer limit used at working point (s,t) 

m(s,t) Multiplier of imaginary traffic used at working point (s,t) 

Advantageously, these state variables vary linearly with the traffic, so a new 
connection can simply add its contributions to them. 

Finding the point (s,t) which optimizes the approximation to the probability of 
25 loss with the actual buffer size limit, service, and traffic, in real time is generally not 
feasible. Instead, the present invention uses a predetermined working point (s,t) and 
adds some imaginary traffic from a design traffic mix to the actual traffic in order to 



loake the working point the optimizing point for the sum of real and imaginary traffic. 
This makes it sufficient to keep track of the state variables rather than of whole 
functions. 

The approximation of the invention uses the state variables to produce the 
following three parameters: (1) an approximation q = z - log(c) to the logarithm of the 
probability of loss; (2) a buffer size limit B; and (3) a multiple m of imaginary traffic 
from a design mix. 

The traffic on all connections is admissable if the four conditions are satisfied. 
The first condition is that q be less than or equal to the QOS log loss requirement. The 
second condition is that B be less than or equal to the limit set by available buffer space 
and QOS delay requirements. The third condition is that m is nonnegative. The fourth 
condition is that the mean input rate of real plus imaginary traffic exceeds the mean 
service rate by no more than admitted by the QOS loss requirement. 

If none of the first three conditions is satisfied, the traffic is inadmissible. If 
some, but not all, of the conditions are satisfied, the algorithm cannot make a 
determination, and therefore, does not admit the traffic. 

To find a good working point, the above calculations are performed off-line for a 
great number of points (s,t), and a working point is picked from a region in the (s,t)- 
plane that performs best with the design traffic mix. If it is desired to widen the 
tolerance against deviations from the design mix, such may be possible by choosing a 
working point that is not optimal for the design traffic mix. If it is desired to handle 
several design traffic mixes, a good working point can be selected for each design mix, 
state variables maintained for each working point, and calculate q, B, and m calculated 
at each working point at connection set-up. If no working point rejects the traffic, and 
at least one working point admits the traffic, the traffic is deemed admissible. 

The present invention applies, e.g., to a single queue and server. Moreover, the 
invention can be generaUzed to multiple queues and servers. 

The present invention has a wide range of applicability, working both on a 
mixture of encoded speech and leaky bucket shaped traffic as well as other types of 
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traffic. The traffic may be of fluid or discrete arrivals type. In examples with known 
exact answer considered so far the method has admitted between 97% and 100% of 
what is theoretically admissible. Although more complex than a pure efficient 
bandwidth method, the invention is implementable in real time. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features, and advantages of the invention will 
be apparent from the following more particular description of preferred embodiments as 
illustrated in the accompanying drawings in which reference characters refer to the 
same parts throughout the various views. The drawings are not necessarily to scale, 
10 emphasis instead being placed upon illustrating the principles of the invention. 

Fig. 1 is a graph illustrating causes of underestimating loss. 

Fig. 2 is a graph illustrating bounding of x - B by ke^^. 

Fig. 3 is a schematic view illustrating a high-priority queue H and a low-priority 
queue L served by a single server. 

15 Fig. 4 is a schematic view illustrating N queues of equal priority served by a 

single server. 

Fig. 5 is a graph showing a simulation of queue length distributions. 

Fig. 6 is a schematic view of an example of queues with shared buffer and 
independent servers. 

20 Fig. 7 is a schematic view of a queueing system in a spatial switch of a wideband 

CDMA telecommunications system. 

Fig. 8A is a graph of a log moment generating function; Fig. 8B is a graph of the 
derivative with respect to s of the log moment generating function of Fig. 8A; Fig. 8C is 
a graph of the derivative with respect to t of the log moment generating function of Fig. 
25 8A. 
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Fig. 9 is a graph showing ON-OFF periodic arrivals. 

Fig. 1 OA is a graph showing the log moment generating function, and Fig. lOB 
and Fig. IOC its derivatives with respect to s and t, respectively, for ON-OFF periodic 
source with R=l, T=l, Ton=0.2. 

5 Fig. 1 1 A is a graph showing the log moment generating function, and Fig. 1 IB 

and Fig. 1 IC are graphs showing the derivatives of the log moment generation function 
of Fig. 1 lA with respect to s and t, respectively, for discrete periodic source with T=l, 

Fig. 12A is a graph showing the log moment generating function, and Fig. 12B 
10 and Fig. 12C are graphs showing the derivatives of the log moment generation function 
of Fig. 12A with respect to s and t, respectively, for a Poisson source with T=l, R=l. 

Fig. 13A is a graph showing the log moment generating function, and Fig. 13B 
and Fig. 13C are graphs showing the derivatives of the log moment generation function 
of Fig. 13A with respect to s and t, respectively, for an ON/OFF Markov fluid source 
15 witii T=l/Qi,2+1/Q2a=l, Ton=l/Qi.2=0.2, R=l. 

Fig. 14A is a graph showing the log moment generating function, and Fig. 14B 
and Fig. 14C are graphs showing the derivatives of the log moment generation function 
of Fig. 14A with respect to s and t, respectively, for always ON or OFF fluid source 
with PoN = 0.2 and R = 1. 

20 Fig. 15 A and Fig. 15B are graphs showing contours of q, B, and m for traffic 

cases 1, 10, respectively. 

Fig. 16A - Fig. 16D are graphs showing contovffs of q, B, and m for traffic cases 
2, 3, 4, 5, respectively. 

Fig. 17 A - Fig. 17B are graphs showing contours of q, B, and m for traffic cases 
25 6, 7, 8, and 9, respectively. 
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Fig. 18A - Fig. 18D are graphs showing contours of q, B, and m for traffic cases 
11, 12, 13, and 14, respectively. 

Fig. 19A - Fig. 19D are graphs showing contours of q, B, and m for traffic cases 
15, 16, 17, and 18, respectively. 

5 Fig. 20 is a graph showing contours of q, B, and m for a 19* one connection 

traffic case. 

Fig. 21 A - Fig. 2 ID are graphs showing contours of q, B, and m for traffic cases 
20, 29, 39, 48, respectively. 

Fig. 22A - Fig. 22D are graphs showing contours of q, B, and m for traffic cases 
10 21, 22, 40, 41, respectively. 

Fig. 23 A - Fig. 23D are graphs showing contours of q, B, and m for traffic cases 
23, 24, 42, 43, respectively. 

Fig. 24A - Fig. 24D are graphs showing contours of q, B, and m for traffic cases 
25, 26, 44, 45, respectively. 

15 Fig. 25A - Fig. 25D are graphs showing contours of q, B, and m for traffic cases 

27, 28, 46, 47, respectively. 

Fig. 26 A - Fig. 26D are graphs showing contours of q, B, and m for traffic cases 
30, 31, 49, 50, respectively. 

Fig. 27A - Fig. 27D are graphs showing contours of q, B, and m for traffic cases 
20 32, 33, 51, 52, respectively. 

Fig. 28A - Fig. 28D are graphs showing contours of q, B, and m for traffic cases 
34, 35, 53, 54, respectively. 

Fig. 29 A - Fig. 29D are graphs showing contours of q, B, and m for traffic cases 
36, 37, 55, 56, respectively. 



Fig. 30- Fig. 30 are graphs showing contours of q, B, and m for traffic cases 36, 
38, 57, respectively. 

Fig. 3 1 is a schematic view of an node according to an embodiment of the 
invention. 

Fig. 32A is a schematic view of an node entity which includes a node main 
processor. 

Fig. 32B is a schematic view of an node entity which serves as an extension 
terminal. 

Fig. 33 is a graph showing a probabiUty density function. - 

DETAILED DESCRIPTION OF THE DRAWINGS 
1 Introduction 

Link Admission Control LAC is part of Connection Admission Control CAC. 
LAC examines whether resources are available for a new connection on a link. 
Resources are channel numbers, buffer space and bandwidth. The present invention 
considers LAC with respect to bandwidth and buffer usage only. 

A link may contain several queue-and-server systems, for instance, egress queues 
in the sending node and ingress queues in the receiving node. LAC must check for 
resources in every such system. The present invention illustrates an exemplary system. 

In the above regard. Fig. 31 shows a representative node 20 with which the 
connection admission control (CAC) technique of the present invention can be 
employed. The particular node 20 which serves as an illustration of the connection 
admission control (CAC) of the present invention is an Asynchronous Transfer Mode 
(ATM) node 20 comprising a switch core 24. Switch core 24 which has plural switch 
core ports, four of the switch core ports being shown as switch core ports 26A - 26D in 
Fig. 31. A node entity 30, also known as a device board, is connected to each of the 
switch core ports. Fig. 31 shows node entity 30A being connected by a bidirectional 
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link 32A to switch core port 26A; node entity 30B being connected by another 
bidirectional link 32B to switch core port 26B; and so forth. It should be understood 
that more than four node entities 30 can be, and typically are, connected to 
corresponding ports 26 of switch core 24, but that only four node entities 30 are shown 
5 for sake of simplification. 

Each node entity 30 performs one or more functions and has, among other 
components hereinafter described, a processor mounted thereon. One of the node 
entities 30, particularly node entity 30A, has a node main processor which generally 
supervises operation of the entire ATM node 20. The other node entities 30, such as 
10 node entities 30B - 30D, have entity processors SOB - SOD, respectively, also known as 
board processors. 

In the particular embodiment shown in Fig. 31, each of node entities 30B - 30D 
serve as extension terminals. Having such function, the node entities 30B - 30D are 
connected by physical Unes or links to other ATM nodes. For example, node entity 30B 
15 is shown as having four physical lines 60B-1 through 60B-4 to other (unillustrated) 

ATM node(s). Although not necessarily labeled, in Fig. 31, the other node entities 30B 
and 30C also have four physical lines extending to other (unillustrated) ATM node(s). 

In general, the ATM node 20 serves to route ATM traffic cells between physical 
lines 60 which connect ATM node 20 to other ATM nodes. For example, ATM traffic 

20 cells incoming to ATM node 20 on physical line 60B-1 can be routed by switch core 24 
to be outgoing from ATM node 20 on physical hne 60C-1. The entity processor of each 
node entity 30 plays a significant role when establishing ATM connections to/from that 
entity. In case of an extension terminal (ET) entity, the estabhshing of an ATM 
connection between a physical line and another node entity (e.g., another extension 

25 terminal or any other type of node entity) is performed by setting up a translation table 
row (hosted in the ATM line module), one for each direction. In the ingress direction, 
the translation assigns an internal VPWCI and an addressee switch port for each 
utilized VPI/VCI on the physical link. The addressee switch port is used to route each 
cell to the right switch port (i.e., node entity). In the egress direction, the translation 

30 assigns the VPWCI to be used on the physical link for each VPWCI used internally 
between two node entities. When actually transmitting cells on the connection, only 
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hardware is involved (e.g., no processors perform any tasks concerning cell transfer). 
In the case of any other type of node entity (e.g., an entity that terminates an ATM 
connection), the principles discussed above apply except for the egress direction in 
which no external VPWCI is assigned. Instead, a termination point (software entity of 
5 the processor) is utilized. 

As mentioned above, the overall operation of ATM node 20 is managed by node 
main processor 40. In order to communicate with the node entities 30, and particularly . 
with the entity processors 50 of the respective node entities 30, certain control paths 
must be estabUshed between node main processor 40 and the entity processors 50 so 

10 that the processors can communicate with one another. The conmiunication is 

performed by cells which are transmitted over the control paths established between 
node main processor 40 and the various entity processors 50. Establishment of these 
control paths is understood with reference to U.S. Patent Application SN 09/249,785 
filed Febmary 16, 1999, "Establishing Internal Control Paths in ATM Node," which is 

15 incorporated herein by reference. 

Fig. 32A shows an example node entity 30A at which node main processor 40 is 
situated. The node entity 30 of Fig. 32A includes a switch port interface module 
(SPIM) 30A-1 which is connected by bidirectional Unk 32A to switch core 24. The 
switch port interface module (SPIM) 30A-1 is connected to bus 30A-2, which is 
20 preferably a UTOPIA standard bus. The node main processor 40 is connected by bus 
30A-2 to switch port interface module (SPIM) 30A-1. As discussed herein, and 
particularly illustrated in Fig. 7, the switch port interface modules (SPIMs) have queues 
which are affected by traffic management. 

Fig. 32B shows an example node entity 30 which serves as an extension 
25 terminal. Like Fig. 32A, the node entity 30 of Fig. 32B has switch port interface 
module (SPIM) 30B-1 and bus 30B-2, with a processor (entity processor 50) being 
connected to bus 30B-2. In addition, bus 30B-2 is connected to ATM line module 30B- 
3. The ATM Une module 30B-3, as hereinafter explained, contains VPWCI translation 
tables used for performing the external/internal VPWCI and internal/external WFWCl 
30 translations described above. The ATM line module 30B-3 is, in turn, connected to line 
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termination module (LTM) 30-4. It is line temiination module (LTM) 30B-4 which is 
connected to the physical lines 60. 

Examples of the components of a node entity 30 are described, for example, in 
the following United States Patent Applications (all of which are incorporated herein by 
5 reference): U.S. Patent Application SN 08/893,507 for "Augmentation of ATM Cell 
With Buffering Data"; U.S. Patent AppUcation SN 08/893,677 for " Buffering of Point- 
to-Point and/or Point-to-Multipoint ATM Cells"; U.S. Patent AppUcation SN 
08/893,479 for "VP/VC Look-Up Function"; U.S. Provisional AppUcation Serial No. 
60/086,619 for "Asynchronous Transfer Mode Switch." The structure and operation of 
10 node 20 is also further understood with reference to U.S. Patent AppUcation SN 

09/188,101, "Asynchronous Transfer Mode Switch," and U.S. Patent AppUcation SN 
09/188,265, "Asynchronous Transfer Mode Switch", both of which are incorporated 
herein by reference. 

In the representative node 20 herein described it is node main processor 40 (see 
Fig. 31) which performs the connection admission control (CAC) procedures of the 
present invention. In this regard, references hereinafter to an algorithm refer to 
connection admission control (CAC) procedures executed by node main processor 40. 
It should be understood that the connection admission control (CAC) procedures of the 
present invention could be performed for a telecommunications node by a processor 
located other than on a device board, or by structures other than a processor (e.g., a 
hardwired circuit). Moreover, although the representative node 20 happens to be an 
ATM-based node, the connection admission control (CAC) techniques of the present 
invention are not confined to ATM or to any particular type of telecommunications type 
(e.g., CDMA). Rather, the connection admission control (CAC) techniques of the 
present invention can be employed in or for any telecommunications node where traffic 
management is an issue. 

2 Approximations to Bandwidth, Queue Length and Loss Ratio 

Consider a general single server queue, A(t) denotes the amount of work which 
arrives in the interval [-t, 0) and S(t) the amount served in the same interval. If more 
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work anived than can be processed, the surplus waits in the queue, if possible. The 
workload process is defined by W(0)=0 and 

W(t) = A(t) - S(t) 

and the queue of unprocessed work at time zero is 

Q = sup W(t) t > 0 

In the foregoing, 



A(t) 


is the amount arriving in [-t,0) 


R 


is the mean arrival rate (of stationary arrival process) 


S(t) 


is the amount that can be served in [-t,0) 


c 


is the mean service rate 


W(t) = A(t) - S(t) 


is the workload in [-t,0) 


B 


is the buffer size 


fw(x;t) 


is the probability density of W(t) 


Ploss 


is the probability of loss 


Mw(s;t) = E{e^^^^^} 


Moment generating function of W(t) 


/zw(s;t) = logMw(s;t) 


is the log moment generating function of W(t) 



If the workload exceeds the buffer size, the excess amount of work, W(t)"B , gets 
lost. The mean amount of work lost in [-t,0) is 

RtP^^ it) > JJ(x - B)f^ (x; t)dx 

There are three causes for the inequality (illustrated in Fig. 1) are that: (1) the 
queue may be non-empty at -t; (2) the queue may be empty in part of [-t,0); and (3)the 
workload may have a maximum within [-t,0). 
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Solving with respect to Pioss(t) gives 



Rt 



Now the integral is upper bounded by an exponential bound on the factor (x-B), 



r 0 X < B 
X - B x>B 



5 where s is a positive parameter and k is a constant chosen such that the curve y = ke^'^ 

touches the straight line y = x - B. This happens with k = - e'^^*^^^ , as illustrated in Fig. 

s 

2. This type of bound is an integrated Chemoff bound. 

The upper bound on the average amount of work lost in [-t,0) becomes 

j^{x - B)f^(x;t)dx < 



_1 ^-(i+jfl) 



re"'f^{x\t)dx < 



s 
s 



10 where Mw(s;t) = E{e^^^^^} is the moment generating function of fw(x;t) and )aw(s;t) = 
logw(s;t) is the log moment generating function of fw(x;t). The bound applies with 
equality if and only if fw(x;t) consists of a single delta-impulse at x=B+l/s, where the 
exponential curve touches the line y = x - B. 



The bound is tightest for the s that minimizes p.w(s;t) - sB - 1 - log s. The s that 
15 gives the tightest bound must satisfy 
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An s satisfying this equation does indeed always give a minimum. To prove this, show 
ttjat the second order derivative is positive, 

1 

5 The proof is outlined using a probabiUty shift argument. The function 

fw(x;t)e""'^^***'^ is a probability density function, because it is non-negative and its 
integral over all x is 1 . A random variable with this probability density has mean 

— is\ t) and variance -r-^ yi^ {s\ t) . Since a variance is non-negative, this concludes 
as as 

the proof. 

10 The lower bound on PiossCO and the upper bound on the integral can be combined 

into the following approximation: 

log^^,(0 - q 
where q = ill^(sU) - Bs - 1 - log{Rst) 

and B = — piy,(s;t) 

as s 

15 for all s,t > 0 

The maximum of q over all time interval lengths t is used as an approximation to 
the probability of loss. Assuming derivatives exist, a necessary condition for maximum 
is dq/dt = 0. Differentiating: 



dq = -sdB + — fiy^is; t) - B - - ds -\' — ^{s\ r) - - 

as s J \o t t ) 



r 



dt 
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Keeping the buffer limit constant, dB = 0, and using the minimizing s, we get the 
necessaiy conditions for maximum 

Using the s and t found above, it is found that 

= - 
dB * 

Thus, -s is the slope of q versus B. 

There is no general guarantee that dq/dt = 0 gives a global maximum. Periodic 
arrival processes give several local maxima and minima. 

A necessary condition that the solution found above is a local maximum rather 
10 than a minimum is d^q/dt^ < 0. 
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This condition has not had any practical significance in the examples tested. The 
solutions it excluded did not admit the traffic anyway for other reasons. 

A necessary condition that a local minimiun is also global is that the queue is 
15 stable, that is, the mean arrival rate is less than or equal to the mean service rate, R < C. 
Unfortunately, this simple condition excludes the results with constant arrival rate and 
constant service rate for all s,t > 0, although these results are exact, see Example 1 in 
Section 5.1 hereof pertaining to constant arrival rates. These results are a buffer size of 
B = 0, and an excess arrival rate of l/(st), which gets lost. To allow for this case, the 
20 weaker condition is used: 
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R - C ^ Re' 



This condition is important in conjunction with periodic arrival processes. They 
give local maxima in every period, even under heavy overload conditions; but the 
condition shows that the maxima can not be global. 

The results are summarized as follows: 





log Pu^ = 9 




where q 


= Hy,is\t) - Bs -1- 


- logiRst) 


and 


B = —fi^{s;t) 
as 


_ 1 

s 


and 


dt 


_ 1 
t 


for s,t > 0 and R — 


C < Re'"" 



(1) 



2.1 Constant Service Rate 

Assume that the arrival and service processes are independent, and that the 
service process has a constant rate, C. This gives 
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^w(s;t) = iiA(s;t) + ^s(-s;t) = kiaCsU) - Cst 



q and its differential become 



q = |^A(s;t) - Bs - Cst - 1 • log(Rst) 



fd 1 
dq = sdB - StdC -hi — fiAs; t) - B - Ct - - 

\ds SJ 



ds + 



(d 



1 



V 



— ^i^(<s\t)-Cs'^-\dt 
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Again, to minimize q for constant B, C, and t, the coefficient of ds must be 0, To 
maximize q for constant B, C, and the minimizing s, the coefficient of dt must be 0. 
The s and t determined like this have the following interpretation: 
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-s is the slope of q versus B for constant C 
-St is the slope of q versus C for constant B 
The result for constant service rate C is 







log Ploss = ^ 


where 


q = PL^ (s; t) - Cst - Bs - I- \og(Rst) 




and 


B = -^fi,is;t)-Ct-- 
as s 




and 




for 


s,t > 0 


and R - C < Re"""" 



2J1 General Service 

Solving equations (1) for (s,t) is difficult. Instead, we shall use (1) with the 
actual arrivals and service replaced by something else. 

The "unused typical arrivals" problem formulation in section 2.2. 1 presents a 
solution in a form suitable for implementation in a real-time CAC algorithm. It' s 
suitability comes from a "single point calculation" property. 

The problem formulations in sections 2.2.2-2.2.4 are special cases and variants, 
of this. Despite the solutions look simpler they are less suitable for implementation in a 
real-time CAC algorithm. This is because they do not have the "single point 
calculation" property, but rather a "leaning banana" trouble. 

2.2.1 Unused Typical Arrivals Formulation 

Replace the actual arrivals A(t) by A(t) + U(t), where U(t) is unused typical 
arrival trsiffic. 
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Generate U(t) as follows: Choose a typical set of connections, and call its mean 
arrival rate Rq and its log moment generating function Mo(s;t). Add independent 
connections in the same proportions as in the typical set, but m times as many, to form 
U(t) with mean arrival rate mRo and log moment generating function ^u(s;t) = m|io(s;t). 
Choose m such that the traffic satisfies equations (1). 

The result for general service S(t) with unused typical arrivals is 



log Pioss = <1 

where q = ^i^{s\t) - tiy^ {s\ t) + n{ ii^{s\ t)- s-^ ix^{s\ t) \ - log((/? + mR^)st) 

as \ as J 

d d \ 

and iB = — /iv|.(^;/) + m — /io(j;0 - - . (3) 

as as s 

d / 1 
— a^{s\t) - - 

and m = -2L£^ ^ 

d t 



for s,t>0 and R + mR^ - C < (/? + mRQ)€ 



fnwx 



It should be noted that m=0 corresponds to the maximum admissible load of real 
traffic. If m is non-negative, the service S(t) is sufficient, otherwise S(t) is insufficient. 
B and q in (3) are not the actual buffer size and log loss probability, but the buffer size 
and log loss probabiUty that would arise if the unused typical arrivals were added to the 
actual arrivals. 

The terms m and B in (3) are linear functions of the derivatives of ^iw(s;t). q is 
the sum of a Unear and a logarithmic terai., 



= z - log c 



d f d ^ 

where z = jU^(^; r) - s—^fiy^(s; t) + rri\^fi^{s\ t) - s—^fi^is; t)j - log(st) 



(4) 



and c = /? + /n/?o 
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A CAC algorithm can keep track of c, z, B, and m by adding contributions from 
new connections and subtracting contributions from connections being cleared. 

The "unused typical arrivals" formulation of the problem has a "single point 
calculation" advantage. If the real arrivals come from a typical traffic mix, then the 

5 mean rate and log moment generating function of A(t)+U(t) do not vary with the real 
traffic load. This means that a point (s,t) that admits the traffic at the maximum 
admissible load also admits at all lower loads. If the real traffic deviates from the 
typical mix the sum of real and imaginary traffic may still be sufficiently similar to the 
typical mix to allow a single point (s,t) to admit over a range of real traffic mixes and 

10 loads. 

2.2.2 Unused Arrival Bandwidth Formulation 

Replace the actual arrivals A(t) by A(t)+Ut, where U is the unused arrival 
bandwidth. 

The result for general service S(t) with unused arrival bandwidth U is 







log Pu^^' q 




where 


q = ^is',t) + Ust - Bs - I - log((R + U)st) 




and 


B = -|-AV(^;0 + ut 

as s 






and 


U = -- —^is;t) - - 
s\at tj 




for 


s,t > 0 


and R + U- C<iR + £/)e'- 



If U is non-negative, S(t) is sufficient, otherwise S(t) is insufficient. B and q in 
Equation Block (5) are not the actual buffer size and log loss probabiUty, but the buffer 
size and log loss probability that would arise if the actual arrivals were increased by the 
unused bandwidth. 



4 



■ 
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The log loss probability q in Equation Block (5) is the sum of a linear and a 
logarithmic term, 

q ^ z- log(R + U) 
where z = fJ^{s;t) Ust - Bs - I - log(jr) 

An admission control algorithm can keep track of B, R, U, and z by adding 
5 contributions from new connections and subtracting contributions from connections 
being cleared. 

A drawback of using Equation Block (5) for real-time admission control is that 
the (s,t)-region admitting the last connection does not always overlap the (s,t)-region 
admitting the first connection. In (s,t,load)-space the region admitting the loads may 
10 look like a 'leaning banana''. Thus, it is not sufficient just to calculate Equation Block 
(5) in a single point (s,t). 

2.2.3 Unused Service Bandwidth Formulation 

Replace the actual service S(t) by S(t)-Ut, where U is the unused service 
bandwidth. 

15 This gives the same workload on the queue as in the unused arrival bandwidth 

formulation. U and B are the same as before, but q is greater when U>0. At full load, 
when U=0, the two unused bandwidth formulations give the same result. 

The result for general service S(t) with unused service bandwidth U is 
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log Pto„ = Q 



where q = + 1/5/ - 5^ - 1 - log(/?5r) 




(6) 



and i7 = — ^{s\t) 

5 V (9f r 




for 



sj>0 



and 



/? " (C - f/) < Re 



If U is non-negative, S(t) is sufficient, otherwise S(t) is insufficient. B and q in 
(6) are not the actual buffer size and log loss probability, but the buffer size and log loss 
probabihty that would arise if the actual service was reduced by the unused bandwidth. 

The set of points (s,t) admitting the arrivals differs between the two unused 
bandwidth formulations. If Equation Block (6) admits the arrivals at (s,t), then so does 
Equation Block (5); but the converse is not true. If Equation Block (5) admits the 
arriv£ds at (s,t), the Equation Block (6) may admit or be uncertain. 

The log loss probability q in Equation Block (6) is the sxmi of a linear and a 
logarithm term, 



where z = Mh^C-^' 0 + Ust - Bs - 1 - log{st) 

B, R, U, and z are all linear functions of the log moment generating function of 
the workload W(t). An admission control algorithm can keep track of B, R, U, and z by 
adding contributions from new connections and subtracting contributions from 
connections being cleared. 

The unused service bandwidth formulation suffers from the same "leaning 
banana" trouble as does the unused arrival bandwidth formulation. 



<7 = z - log /? 
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2.2.4 Arrivals Multiplier Formulation 

Imagine increasing the arrivals k times by adding independent connections from 
the same traffic mix as the actual connections. The resulting arrivals have mean rate 
and log moment generating function 

^ kR 

Select the multiplier k such that the workload resulting from imaginary arrivals 
and actual service satisfies Equation Block (1). The result is 



log = 9 

where q = kfi^^s; t) -h fi^{-s;t) - Bs - I - log(kRst) 



and 



as as s 

and A:=--^ ^ 



dt 



for 



/ > 0 



kR - C <: JtRe^™* 



(7) 



If k^l, B^Bo, q^qo at some (s,t), the connections are admissible. If k<l, B>B^, 
10 q>qo at some (s,t), the coimections are inadmissible. 

The admission control algorithm of the present invention keeps track of four 

3 d 

state variables R, /i^ {s; r), -r- fx^ (s\ r), — /i^ (5; t) by adding contributions from new 

as at 

connections. The algorithm can calculate q, B, and k from these state variables. 

The arrivals multiplier formulation suffers from the same problem as the unused 
15 bandwidth formulation, that the point (s,t) that admits the last coimection at maximum 
load will in general not admit the first connection. 
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2.3 Asymptotic Expressions for Large t 

Assume that the workload process W is stationary and ergodic. Assume that it 
satisfies the conditions of the Gartner-Ellis theorem; that is, assume that the asymptotic 
log moment generating function 

exists and is finite for all real s, and that hw is differentiable. That hw is convex, 
positive and increasing for s>0 can be direcdy verified. This means that for large t the 
following approximation can be used: 

For constant rate and Poisson process, this product form of KLw(s;t) holds for all t. 
For other interesting processes, e.g. leaky bucket shaped traffic, it is poor for small t. If 
|iA(s:t)=hA(s)t and the service process has constant rate C, the results becomes 



log - q 

where q = h^is)t - Cst - Bs - 1 - logiRst) 

and 5 = h^' is)t - a - i 

s 



and 



S 



r 



1 



for 



> 0 



i? < c + — 

St 



3 Delay 

3.1 Constant Service Rate 

For a given pair (s,t), Equation Block 3 gives a buffer size B and the logarithm q 
of the probabiUty of exceeding B. Data arriving to a B long queue served at constant 
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rate C will leave the queue after a delay of D=B/C. If the imaginary traffic is positive, 
the delay of the real traffic will be smaller. Hence, the (1-e^) fractile of the delay is at 
most B/C, and 



P{D >^)^e' 



(8) 



where the equality sign applies when the system is fully loaded. 



4 Multiple Queues 

4.1 Multiple Priorities, Single Server 



Consider a system consisting of two queues, H and L, served by a single server 
as shown in Fig. 3. The server serves queue H whenever that queue is not empty. The 
10 server serves queue L when queue H is empty and queue L is not. Queue H sees the full 
service of the server. Queues H and L together also see the full service of the server. 
Queue L sees what queue H leaves unused. The log moment generating function of the 
service offered to queue L is determined as follows. 

The workload on queue H is WaCt) = AH(t)-S(t). The service offered to queue L 

15 is 

S,(r) = (S(r) - A^iOr = "Wnit) 

\W^(0 for W„(r)<0 

where W H{t)\ 

[ 0 for W^(t)>0 

ft 

From this, the log moment generating function of the service offered to queue L 
is determined: 



« 



r 
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This is a function of /x^^ (s; t) , but unfortunately, a non-linear one. This means that 

contributions can not just be added to it from every high-priority connection. A 
practical approximation is to replace Wh' by Wh and use 

fis, (5; 0 - /V^ i's; t) (9) 

5 This is equivalent to loading queue L by the total traffic A and offering it the full 

service S. It is a good approximation if the probability of positive workload on queue H 
is small. 

4.1.1 Admission Control for Multiple Priorities^ Single Server, Shared 
Buffer 

to 

10 In this scenario, queues H and L share a common buffer of limited size. 

Equation Block (3) is used with workload Wh = Ar - S to evaluate qH at a pair (s,t). 

If qH is admissible, use Equation Block (3) with workload W = A - S. Evaluate 
the overall variables m, B, and q at (s,t). Find the low-priority log loss probability as 
the solution of the following loss rate equation. 

15 Re^ = Rf^e"^" + R^e""' 

Re^— R e^" 

If qL'^O, the solution is qL = log " . To prevent division by 0, check 

that 

Re' - R^e""" < R^e'*'-^' (10) 

instead of checking qL^qimax directly. 

20 Show that all of m, B, qn, and qL are admissible or inadmissible at (s,t), if 

possible. 
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There may be more than two priority levels. Evaluate of the each lower level 
as above using the q of the total traffic on this and all higher levels. At the lowest level, 
show that m, B, Qh, and all are adnussible or inadmissible at (s,t), if possible. 

4.1.2 Admission Control for Multiple Priorities, Single Server, Individual 
5 Buffers 

In this scenario, queue H has a buffer of limited size, and queue L has a buffer of 
limited size. Equation Block (3) is used with workload Wh=Ah-S. Show that all of the 
high priority variables mH, Bh, and qn are admissible or inadmissible at a pair (sH,tH), if 
possible. 

10 If mH, Bh, and qn are admissible, use Equation Block (3) with workload Wl=Al- 

Sl, approximating Sl as in Equation Block (6). Show that all of the low-priority 
variables mL, Bl, and qL are admissible or inadmissible at a pair (SL,tL), if possible. 

There may be more than two priority levels. Repeat the evaluation for layer L 
for each lower priority level using the sum of traffic on all higher layers for Ah- 

15 4.2 Single Priority 

Consider a system consisting of N queues served by a single service as illustrated 
in Fig. 4. The server serves the non-empty queues according to a schedule giving them 
an equal share of the service. This schedule could be round robin or according to a 
snapshot of queue states taken when a previous snapshot has been processed. The N 
20 queues together see the full service of the server. An individual queue sees what the 
other queues leave unused. The situation looks like the one for a low-priority queue 
seeing the service the high-priority queues leave unused, and we shall use the same 
approximation (Equation 9), restated here for convenience. 

Ms, 0 ^ ti^v„ (-^y; 0 (9) 



25 Here Sl is the service seen by an individual queue, and Wh is the total workload on the 
other queues. Again, this approximation gives the same distribution of the length of one 
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queue as the distribution of the sum of the lengths of all the queues. Fig. S shows a 
simulation result justifying such approximation. 

Fig. 5 shows distributions of queue lengths. The system comprised 21 queues of 
equal, low, priority served by a single server. Equal traffic load on all queues with 
exponentially distributed time between bursts, geometrically distributed number of 
arrivals in a burst, and constant inter-arrival time within a burst. Scheduhng according 
to a snapshot of queue states, queue 0 served first and queue 20 last. In Fig. 5, the 
lower set of curves shows the length of queue 0, which was served first in the snapshot; 
the middle set of ciu^es shows the length of queue 20, which was served last in the 
snapshot; and the top set of curves shows the sxmi of the lengths of all 21 queues. 

4.2.1 Admission control for Single Priority, Single Serl^er, Shared Buffer 

In this scenario, the N queues share a conmion buffer of limited size. This case 
is equivalent to the single queue case. Equation Block (3) is used with workload W=A- 
S. Evaluate the overall variables m, B, and q at a pair (s,t). Show that all of U, B and q 
15 are admissible or inadmissible at (s,t), if possible. 

4.2.2 Admission Control for Single priority. Single Server, Individual 
Buffers 

In this scenario each of the N queues has a buffer of limited size. Let Al be the 
arrivals to the queue considered, and let Ah be the sum of arrivals to the other queues. 
20 Equation Block (3) is used with service Sl=S-Ah. Show that the variables mt, B^, and 
qL of the queue considered are all admissible or inadmissible at a pair (sL,tL), if possible. 

4.3 Queues with Individual Service and Shared Buffer 

The sum of workloads on queues with independent servers does not give a 
correct estimate of the sum of queue sizes. Fig, 6 shows an example illustrating this 
25 point. In Fig. 6, queue 2 is unstable and grows beyond all limits, while the sum of 
workloads on the queues is less than the sum of services. It is postulated that the total 
buffer needed for queues with independent servers is approximately equal to the need of 
the longest queue. A motivation is that the queues are independent and a long queue 



5 
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has a small probability, so when one queue is long, the others are likely to be around 
their mean length which is much smaller. 

4.4 Queues in the SPAS Switch 

» 

Fig. 7 shows a simplified diagram of the queueing system in the spatial switch 
5 SPAS (core 24 in Fig. 31) in a wideband CDMA telecommunications system. As 
mentioned above, details of the spatial switch SPAS beyond those present here are 
discussed in U.S. Patent Application SN 09/188,101, "Asynchronous Transfer Mode 
Switch," and U.S. Patent Application SN 09/188,265, "Asynchronous Transfer Mode 
Switch", both of which are incorporated herein by reference. 

10 The switch core 24 comprises rows for incoming data and columns for outgoing 

data. In the crosspoints between rows and columns there are small buffers which are 
omitted in Fig. 7. The space switch interface modules SPIM contain ingress and egress 
queues towards and from the switch core, respectively. 

The ingress queues in a SPIM are organized as a queue per egress SPIM and 
15 ingress priority. The idea is to store incoming cells in ingress queues and transfer them 
in an ingress priority order towards the switch core when the corresponding crosspoint 
is free. Ingress queue (i, j, p) is located in SPIM i and carries traffic towards SPIM j on 
priority level p. 

The ingress queues in a SPIM share a common buffer memory. Ingress queues 
20 in different SPIMs do not share a buffer. 

The reading of cells from the ingress queues follows the following rules: 

• Cells with the highest priority will be fetched until the queues are empty 
or blocked by occupied crosspoints. 
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Lower priority cells will be read as soon as there are no cells in queues 
with higher priority or that all of the higher priority cells are blocked by 
occupied crosspoints. 
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A round-robin mechanism is used to give equal bandwidth to the different ingress 
queues within the same priority level. 

The switch core has two priority levels. To get a high priority cell through the 
s>yitch core when a low priority cell blocks the crosspoint, a command is sent as a 
special management cell called Plus Priority cell. This cell is terminated in the core and 
any cell in the cross-point buffer gets its priority increased. In this way, the crosspoint 
will be emptied within a predetermined time. 

The egress SPIM serves the crosspoints in its column according to a snapshot 
mechanism. When the SPIM has processed a previous snapshot it takes a new snapshot 
of the crosspoint buffer states. The SPIM empties the non-empty buffers in the 
snapshot in increasing number order. 

The egress queues are organized as a queue per egress priority and outgoing 
branch from the SPIM. The idea is to store incoming cells in queues and transfer them 
in egress priority order. The egress queues in a SPIM share a conmion buffer memory. 
This memory is separate from the ingress buffer memory. 

4.4.1 Ingress Admission Control 

This section describes how to use the "unused typical arrivals" results for 
admission control of the ingress queues. 

Row i and colunm j in the switch core both limit the service of the ingress queues 
in SPIM i towards SPIM j. The ingress admission control algorithm checks that both of 
these potential bottlenecks give sufficient service. 

The ingress queues served by a row share a buffer in a SPIM. Hence it suffices 
to keep track of the total workload per priority level when checking for row service. 

The ingress queues served by a colunm reside in different SPIMs, and therefore 
do not share a buffer. Hence one must keep track on workload on individual queues 
when checking for column service. Queues sharing a buffer in a SPIM do not share 
colunm service. 
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The ingress admission control algorithm for row service uses the following state 
variables for some set of points (s;t), 

CR(p,i) Average rate of real + imaginary traffic on priority levels p 

and higher from SPIM i 

ZR(p,i,s;t) Linear term in long loss probability of traffic on priority 

levels p and higher from SPIM i 

BR(p,i,s;t) Buffer size of traffic on priority levels p and higher from 

SPIMi 

mR(p,i,s;t) Multiplier of unused typical traffic on-priority levels p and 

higher on row i 

The ingress admission control algorithm for column service uses the following 
state variables for some set of points (s;t), 

Cc(p4d) Average rate of real + imaginary traffic on priority levels p 

and higher from SPIM i towards SPIM j 

zc(p4d>s;t) Linear term in log loss probability of traffic on priority 

levels p and higher from SPIM i towards SPIM j 

Bc(p,iJ,s;t) Buffer size of traffic on priority levels p and higher from 

SPIMi towards SPIM j 

nic(pj,s;t) Multiplier of unused typical traffic on priority levels p and 

higher on column j 

The typical traffic has the following characteristics for some set of points (s;t), 

Ro(p,i) Average rate of typical traffic on priority levels p and higher 

from SPIM i 



30 



^o(p4,s;t) Log moment generating function of typical traffic on 

priority levels p and higher from SPIM i 

— fj^ (/?, /, s; t) Partial derivative with respect to s, tj^ical traffic on priority 

a s 

levels p and higher from SPIM i 



— fj^ ip, f , s\ t) Partial derivative with respect to t, typical traffic on priority 

o t 

levels p and higher from SPIM i 

For shortness, the coefficients of are introduced: 

dt 



^0 

ag{p,i,s\t) = in the expression for R 

a^ip, i, s; t) = in the expression for z 

d 

Ugip.i, y, r) = -3^^ in the expression for B 

at 

1 

a„(/7, i, s; t) = -3 in the expression for m 



The initial values for row service, with no traffic through the queues, are 
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ip, i, s; 0 = fl/i ip, i. t) 



/, ^; r) = -C,/ + a^C/^. ^ 0 



r 



(11) 



The initial values for column service, with no traffic through the queues, are 



Cc(P. ^ 7. -Si t) = flg (p, i, s; t) 



f 1"! 

C,5 + - 



( n 

Zcip> ^ 7. 0 = «z(P' ^ -^^ ^y"^ + - ~ logC-^O 



5c(P' ^ J' 0 = -Cjt + (p, i, r, r) 



tJ S 



m,ip,i,J,s,t) = aSpJ,s;t)\C,s^\ 



(12) 



Input to the ingress admission control algorithm about a new connection is 



Priority level 



Ingress SPIM 



J 



Egress SPIM 



Average rate 
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Kia(s;t) 



Log moment generating function of arrivals at a set 
of points (s;t) 
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Partial derivative with respect to s 



Partial derivative with respect to t 



The new connection adds the following contributions to the row state variables. 



Ac^ (p, I, ^; r) = r - (p, /, s\ t) — //^ {s\ t) 

d d 

t^R {p. U s\ t) = Mo U 0 - ^ IT (•^' 0 - ip. I, ^, 0 3" /i« (^; 0 

^ d 
(p. s;t) = — tx^ {s\ t) - (p, s\ t) — u„ (s; t) . 

OS o t 



Am. (p, I, s; r) = -a^{p, f, r) — (y, 0 



for p> Pa. all (r, r) 



(13) 



5 The new connection adds the same contributions to the affected column state 

variables. 



Acc(p, k, J, s; t) ^ r - a^(p, s; t) '^^t^a(^'^ 0 

d d 
as o t 

d d 

t^c^p.kj,s\t) = -^fi^is; t) - a^{p,Us\t)— jL^{s\t) 

OS o t 

d 

^c(P^ 7' 0 = -flm(P' ^'^ 0 -J-^f^a 0 

for p > p^, all all (5-; f) 



(14) 



4.4.2 Egress Queues 
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An outgoing link serves its egress queues at a constant rate. A problem is the 
traffic arriving to the egress queues. This traffic has changed by passing through the 
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ingress queues and the switch core. As an approximation, assume that the traffic has 
not changed. This situation is equivalent to handling ingress queues with row service 
only, see above. 

4.5 The Admission Control Algorithm 

The algorithm initializes the state variables c, z, B, and m, as described in 
Equation Blocks (1 1)-(12). In order to admit a new connection, the algorithm checks 
every affected priority level. If all priority levels admit the connection, it is admissible, 
otherwise not. 

In order to check a queue, the algorithm checks it over one or more time intervals 
t. If no t rejects the connection, and some t admits it, the queue admits the connection 
over the time interval t, otherwise the queue rejects the connection. 

In order to check a queue over a time interval t, the algorithm checks it over one 
or more values of s for the given t. If some point (s;t) rejects the connection, the time 
interval rejects it. If no point (s;t) rejects, and some point (s;t) adepts the connection, 
15 the time interval t admits it. 

In order to check a queue at a point (s;t) the algorithm increments the state 
variables c, z, B, and m, as described in equations (13)-(14). It checks the rate 
condition last in equations (3). If the rate condition is OK, the algorithm calculates q as 
described in equation (4). It corrects the loss probability for losses on higher priority 
20 levels at the same point as described in equation (10). If the corrected loss probabiUty is 
sufficiently small, and the buffer B is sufficientiy small, and the traffic multiplier m is 
non-negative, and the rate condition is OK, then the queue admits the connection at 
point (s;t). If the corrected loss probabiUty is negative, then the queue rejects the 
connection at point (s;t). Otherwise the queue is indecisive at point (s;t). 

25 If the algorithm admits a new connection, it updates all affected c, z, B, and m. 

If the algoritiim rejects a new connection, it reverts to the previous values of all affected 
c, z, B, and m. 
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In order to release an established connection, the algorithm decrements all 
affected c, z, B, and m by the amounts it added when it adniitted the connection. 

5. Examples of Log Moment Generating Fwictions 

5.1 Example 1. Constant Arrival Rates 



The workload in an interval of length t is 



A(t) = Rt 

where R is the arrival rate. The log moment generating function (illustrated in Fig, 8A, 
for constant rate source with R=l) and its derivative with respect to s (illustrated in Fig. 
8B) and its derivative with respect to t (illustrated in Fig. 8C) are 



H^{s;t) = log J%"5(jc - Rt)dx = sRt 



^li^{s\t) = Rt 
as 

d 

— i"^(^;0 = sR 

at 

d ^ 

— ju^is-t) = 0 
as 

d ^ 
asat 



dt 



^^As;t) = 0 



In this case, the asymptotic approximation is exact: 

HA(s;t) = hA(s)t 
where hA(s) = sR 



The approximation with constant service rate C become 
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toss 



where 



9 = 



- logiRst) 



and 



5 = 0 



and 



C 




for 



sj>0 



and 



— < Re^"« 

St 



This is the case where the bound applies with equality, since fw(x;t) consists of a 
single delta-impulse at x = (R"C)t = B+l/s, where the exponential curve touches the Une 
y = X - B. In this case, the approximation is exact. 

5.2 Example 2. ON-OFF Periodic Fluid Source 

The source is ON for Ton and OFF for Toff. The period is Ton + Toff = T. In the 



ON state, the source generates data at the peak rate R — , and in the OFF state, it 



generates no data. This is an extreme behavior acceptable by a leaky bucket regulator 
with mean rate limit R, and bucket size RToff. The phase of this periodic pattern is 
uniformly distributed in [0,T). 

Fig. 9 is a graph showing ON-OFF periodic arrivals. In Fig. 9, A(T,r+t) is the 
amount arriving in [x,t+x). The phase x is uniformly distributed in [0, Ton + Toff). From 
Fig. 9, it can be seen that the density function of A is the sum of two delta-impulses and 
a uniform distribution between them. 



T 



on 



36 

Six - nRT) for t = nT 
aS(,x - X,) + bU(Xj jc, , jCj) + cS{x - Xj) otherwise 

where 
1 



for X, < X < X- 
0 elsewhere 



For t equal to an integer nunfiber of periods, this degenerates to a single delta- 
impulse. For t equal to an integer plus one half number of periods and Ton = Toff, this 
degenerates to a imiform distribution. Table 1 lists the parameters and their partial 
derivatives for 0<t<T. 

TABLE 1 

ON-OFF PERIODIC ARRIVALS, PARAMETER VALUES IN DENSITY 

FUNCTION 





0<t^Ton,Toff 




Toff^t^Ton 


Ton,Toff^t<:T 


a 




off 


off 


* off 
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T 
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0 




on 
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2t 


2T 
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2(T - t) 
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T 




R — t 
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on 


RT 


T 

R — t 
T 

on 


RT 


C 


on 


t - T 

on 


T„ - t 

an 


t - 

on 
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T 


T 


da 


- 1 


- 1 


1 


1 


Tt 


T 


T 


T 


T 



37 







T <t<T a- 


T rc<f<T 

A off ^ "-^ on 




dt 


0 


0 


rX 

T 

on 


R^ 


db 


2 


0 


0 


-2 


dt 


T 






i 


dt 


rI- 

T 

on 


0 


rX 

T 

on 


0 


dc 


- 1 


1 


- 1 


1 


dt 


T 


T 


T 


T 



For t = nT + aT, 0<a<l use 

fA(x;nT+aT) = fA(x-nRT;aT) 

The following auxiliary variables are now introduced: 

y = s(x2-xi) 
z = e"'' 
1-z 

V = 

y 

w = az + bv + c 

with partial derivatives 



» 



and 



dv 
ds 



— * 



dv 

dt 
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iy- 

OS 

d s d s 

d sy yl ds 
Idy 

2 ds 



for y ^ 0 



for y = 0 



^ w d z d V 
fl^ J J d s 



d y ( dxj dx^ 
d \ dt d t 



J 



dz_ _ d_y 
dt d t 



dzl 

d ty 



(l-z) dy 
dt 



y2 

Idy 
2 d t 



for 3? 5t 0 



for y = 0 



dw da dz db ^dv dc 

= Z Cl J- V + 1 

d t dt d t d t dtdt 



The log moment generating function and its derivatives can now be written as 



fiAs;nT -h cO^ = s(x^ + nRT) + log w 



— fi^(s\ nT + aT) = jc^ + nRT + — 
ds 



1 d w 
w ds 



— pi^{s\nT OCT) ^ s 
o t 



dxn I dw 



dt 



+ — 
w 



dt 



Fig. 1 OA is a graph showing the Log moment generating function, and Fig. lOB 
and Fig. IOC its derivatives with respect to s and t, respectively, for ON-OFF periodic 
source with R=l, T=l, Ton=0.2. 
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5.3 Example 3. Discrete Periodic Arrivals 

An amount of RT arrives at times ^ = nT + aT, where n is integer and ae [0,1) is 
constant. The phase aT is uniformly distributed in [0,T). 

In an interval of length r = nT + aT arrives an amount A(t) of either (n+l)RT or 

nRT, 

P(A(nT + aT) = nRT) = 1-a 
P(A(nT + aT) = (n+l)RT) = a 

The log moment generating function, its derivatives, and the asymptotic log 
moment generating functions are 

fx^(s; nT + oO^ = sR(n + i)T + log((l - a)e'"^ + a 



f^As;nT + oT) = Rin + 1)T - 



(1 - a)e~"^ + a 



r) 1 1 — 



HAs\nT-¥Cir) = 



ds^^ ((1 - a)e-"^ + a) 

(j; nT + Cu) — 



dtds HI - a)e-"^ + a) 

H. {s\ nT + aT) = - ^ ' 



dt^ T^{(l-a)e-"^ + a) 

h. ( j) = sR 



Notice that the asymptotic log moment generating function is the same as for a 
constant rate process. Using this, the asymptotic log moment generating function in the 
approximations is thus equivalent to using a constant rate model of the arrival process. 
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Fig. 1 1 A is a graph showing the log moment generating function, and Fig. 1 IB 
and Fig. 1 IC are graphs showing the derivatives of the log moment generation function 
of Fig. 11 A with respect to s and t, respectively, for discrete periodic source with T=l, 
R=l. 

5.4 Example 4. Poisson Arrivals 

The mean interval time is T. Each arrival contributes an amomit of RT. The 
probability of k arrivals in time t is (e"^(t/T)'')/k!. The log moment generating function 
and its derivatives of the arrival process are (as explained in A. Papoulis, "Probability, 
Random Values, and Stochastic Processes", McGraw-Hill, 1965), 



d 
ds 



— Ma (-y; 0 = Rte"^ 



d 1 .«T 



ds 

£_ 

dsdt 
d^ 



2 f^A^S't) = R'Tte 



= 0 

Also in this case, the asymptotic approximation is exact: 

tXA(s;t) = hA(s)t 
where hA(s) = (e'^'^-l)/T 



Fig. 12A is a graph showing the log moment generating function, and Fig. 12B 
and Fig. 12C are graphs showing the derivatives of the log moment generation function 
of Fig. 12A with respect to s and t, respectively, for a Poisson source with T=l, R=l . 
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5.5 Examples. Markov Fluid Arrival Process 

A source is called a Markov fluid if its time-derivadve is a function of a 
continuous-time Markov chain on a finite state space. Let 1, . . m be the state space 
and 

Qa = [Qij] th^ irreducible transition rate matrix 

Qij the transition rate from state i to another state j 

qii = - Qij minus the sum of transition rates from state i 

n = (tci, . . Ttm) the row vector of steady-state probabilities, is the 

solution of kQa = 0 and Til = 1 

Ra = diag(Ru ...» Rm) diagonal matrix of arrival rates 

Rj arrival rate in state i 

1 a column of I's 

X(X) The largest real eigenvalue of a matrix X 

The log moment generating function (see Kesidis et al, "Effective Bandwidths 
for Multiclass Markov Fluids and Other ATM Sources", IEEE Trans. Networking, Vol. 
1, No. 4, pp. 424-428, Aug. 1993), and its partial derivatives and asymptotic log 
moment generating function are 



10 
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— a (s- = ^ A— e<c* II 

The partial derivative with respect to s is not simple in general. To find it, 
diagonalize 

Qa+sRa = V(s)D(s)V-Hs) 
where D(s) is a diagonal matrix of eigenvalues of Qa+sRa 
and V(s) is a matrix of column eigenvectors 

Now use 



CO 



and find after some calculation 



ds ds ds 



5.5.1 HIGH/LOW Rate Source 



A HIGH/LOW rate source is a special case with two states, 
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Qa = 



[22.. -(22.. J 
/?, 0 



(transitions / s) 



0 R. 



(bit / s) 



and an ON/OFF source is special case of this with R2=0. 



The steady state probabilities of the source are 



n = 



Qia +Q2.1 e,.2 +(22.1. 



Direct calculation gives the eigenvalues 



h 



(j) = A, = a + Va^ - 4b^ 
Aj = - a - Va^ - 46 j 



where a = G, j + gj, - ^(i?^ + R^) 
and 6 = s^R^R^ - siQ^^^R, + gj../?,) 



Diagonal eigenvalues and column eigenvectors matrices are 



Vis) = 



D(s) = 



\ 0 



L 



Qi.2 -sRt+\, -sR.+A^ 



Derivatives with respect to s are 
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dP 

ds 



d\ 0 

ds dXj, 
0 lb 



dV 
ds 



0 



-R, + 



0 



d\ 
ds 



ds ^ 



where 



dX, 
ds 



da 
ds 



r 



V 



da 
ds 



-2 



ds 



\ 



J 



2Va^ - Ab 



and 



ds 



da 
ds 



^ da ^db^ 
a — -2 



ds 



ds) 



and where 



da 
Is 



- - Ri 



and 



db 
ds 



= 25/?,/?2 -Qu^2 'Q2aRi 



Fig. 13A is a graph showing the log moment generating function, and Fig. 13B 
and Fig. 13C are graphs showing the derivatives of the log moment generation function 
of Fig. 13A with respect to s and t, respectively, for an ON/OFF Markov fluid source 
with T=l/Qi,2+1/Q2,i=l, Ton=l/Qi.2=0.2, R=L 



5.6 Always ON or OFF Fluid Arrival Process 



This is a limiting case of a process with very heavy tails on the holding time 
distributions. The source is ON always with probability Pon» generating data at rate 
R/Pqn, or it is OFF always with probability Poff=1-Pon» generating no data. The 
probability density function of the amount of data received in time t is 
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The log moment generating function and its derivatives are 



tJL^ {s\ t) = log 



s 1 

P P e 



J 



ds 



Rt 



p e P 

*^OFF^ ^ ^ON 



dt 



sR 



R 

-J 1 

^OFF^ ^ ''ON 



sRt 



OAT 



flew 



P e 

* OFF^ 



+ P. 



ON 



J 



sRt 



sRt 



P^..e ''^ + + sRt e 



OFF' 



ON 



dsdt 



fiAs\t) = R 



ON 



/ sRt \ 

flow 



Pqff^ 



+ P, 



ON 



J 



sRt 



ON 



dt 



ON 



f sRt^ \ 

Pqn 



P e 

' OFF^ 



+ P 



ON 



J 



Fig. 14A is a graph showing the log moment generating function, and Fig. 14B 
and Fig. 14C are graphs showing the derivatives of the log moment generation function 
of Fig. 14A with respect to s and t, respectively, for always ON or OFF fluid source 
with PoN = 0.2 and R = 1. 
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6 Examples of Loss, Buffer, and Bandwidth Consumption 

This section defines some traffic classes spanning a wide range of bandwidth and 
burstiness. It uses traffic cases spanning a wide range of link loads to illustrate regions 
of (s,t) adniitting or rejecting a traffic case. Table 2 lists traffic classes. 

5 Fig. 15 A -Fig. 30B show contours of q = In 10'^ = -18.4, B =3 Mbit, and m = 0 

bit/s, and m = mmax, for differing traffic cases and the region(s) bounded by these 
contours admitting or rejecting the particular traffic case. Fig. 15A -Fig. SOB assume a 
link bandwidth of Co = 150 Mbit/s. The typical traffic used is the actual traffic itself 
unless otherwise stated. 



Table 2. Traffic Classes 


Class 
ID 


Source model 


Average rate, 
R (bit/s) 


Average 
period, T(s) 


Probability of 
ON, PoN 


1 


Constant rate 


10* 






2 


ON-OFF periodic 


10* 


0,01 


0,01 


3 


Discrete periodic 


10^ 


0,01 




4 


ON-OFF Markov 


10* 


0,01 


0,01 


5 


Poisson 


10* 


0,01 




6 


ON-OFF periodic 


10* . 


100 


0,01 


7 


Discrete periodic 


10^ 


100 




8 


ON-OFF Markov 


10* 


100 


0,01 


9 


Poisson 


10* 


100 




10 


Always ON or OFF 


10* 




0,01 


11 


ON-OFF periodic 


10' 


0,01 


0,1 


12 


Discrete periodic 


10' 


0,01 




13 


ON-OFF Markov 


10' 


0,01 


0,1 


14 


Poisson 


10' 


0,01 
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15 


ON-OFF periodic 


10' 


100 


0,1 


16 


Discrete Periodic 


10' 


100 




17 


ON-OFF Markov 


10' 


100 


0.1 


18 


Poisson 


10' 


100 




19 


Always ON or OFF 


10' 




0,1 



4 



48 

Table 3 lists numbers traffic cases with one connection. 



Table 3 


Case ID 


Description 


1 


1 Class 1 


2 


1 Class 2 


3 


1 Class 3 


4 


1 Class 4 


5 


1 Class 5 


6 


1 Class 6 


7 


1 Class 7 


8 


1 Class 8 


9 


1 Class 9 


10 


1 Class 10 


11 


1 Class 1 1 


12 


1 Class 12 


13 


1 Class 13 


14 


1 Class 14 


15 


1 Class 15 


16 


1 Class 16 


17 


1 Class 17 


18 


1 Class 18 


19 


1 Class 19 



Fig. 15 A illustrates the foregoing one connection traffic case 1; Fig. 15B 
illustrates one connection traffic case 10. Both Fig. 15A and Fig. 15B show contours of 
5 q, B, and m, for one connection, average rate lOkbit/s. Fig. 15A shows case 1 (constant 
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rate); Fig. 15B shows case 10 (always ON or OFF, Pon=0.01). No B>3Mbit is found 
for m<mmax; no m<0 is found for m<mmax- 

Fig. 16A - Fig. 16D illustrate traffic cases 2, 3, 4, 5, (showing contours of q, B, 
5 and m) for one connection, average rate lOkbit/s, average period 10 ms. Fig. 16A 

shows case 2 (ON-OFF periodic, Pon=0.01); Fig. 16B shows case 3 (discrete periodic); 
Fig. 16C shows case 4 (ON-OFF Markov, Pon=0.01); and Fig. 16D shows case 5 
(Poisson). 



Fig. 17A - Fig. 17B illustrate traffic cases 6, 7, 8, 9 (showing contours of q, B, 
10 and m) for one connection, average rate lOkbit/s, average period 10 s. Fig. 17A shows 
case 6 (ON-OFF periodic, Pon=0.01); Fig. 17B shows case 7 (discrete periodic); Fig. 
17C shows case 8 (ON-OFF Markov, Pon=0.01); and Fig. 17D shows case 9 (Poisson). 



Fig. 18A - Fig. 18D illustrate traffic cases 11, 12, 13, 14 (showing contours of q, 
B, and m) for one connection, average rate lOkbit/s, average period 10 ms.. Fig. 18A 
15 shows case 1 1 (ON-OFF periodic, Pon=0.01); Fig. 18B shows case 12 (discrete 

periodic); Fig. 18C shows case 13 (ON-OFF Markov, Pon=0.01); Fig. 18D shows case 
14 ( Poisson). 



Fig. 19A - Fig. 19D illustrate traffic cases 15, 16, 17, 18 (showing contours of q, 
B, and m) one connection, average rate lOkbit/s, average period 10 s. Fig. 19A shows 

20 case 15 (ON-OFF periodic, Pon=0.01); Fig. 19B shows case 16 (discrete periodic); Fig. 
19C shows case 17 ( ON-OFF Markov, Pon=0.01); Fig. 19D shows case 18 (Poisson). 
Traffic case 16 of Fig. 19B is clearly inadnaissible. A burst of 1 Gbit arrives 
instantaneously to a buffer of size 3Mbit, which looses 99.7% of the arriving data. Yet 
Equation Block (3) admits the connection at some points (s,t) with t>l period = lOO s. 

25 The reason for this is that the workload over long periods looks acceptable, while it is 
unacceptable over short periods (see also Fig. 1). The global maximum of q occurs in 
the first period of t. The admitting points found are at local maxima in later periods and 
thus false solutions. The conclusion of this is that one must choose the working point 
(s,t) for periodic traffic within the first period. Another peculiarity of traffic case 16, is 

30 that no working point (s,t) is able to reject the case. In the region below the contour 
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labeled "-18.4*' in Fig. 19B, q and B are far from admissible, while m is close to 0; but 
m>0 everywhere. 

Fig. 20 illustrates traffic case 19 (showing contours of q, B, and m) for one 
connection, average rate lOMbit/s. Case 19 is always ON or OFF, Pon=01. 

5 Table 4 lists numbers traffic cases (cases 20 - 38) with the maximum number of 

connections from a single traffic class admitted by Equation Block (3). 



Table 4 


Case ID 


Description 


20 


15000 Class 1 


21 


14999 Class 2 


22 


14999 Class 3 


23 


14996 Class 4 


24 


14998 Class 5 


25 


10409 Class 6 


26 


84 Class 7 


27 


10342 Class 8 


28 


81 Class 9 


29 


9508 Class 10 


30 


14 Class 1 1 


31 


14 Class 12 


32 


9 Class 13 


33 


1 1 Class 14 


34 


1 Class 15 


35 


0 Class 16 
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Table 4 


Case ID 


Description 


36 


1 Class 17 


37 


0 Class 18 


38 


1 Class 19 



Table 5 lists numbers traffic cases (cases 39 - 57) with the minimum number 
connections from a single traffic class not admitted by Equation Block (3)* 



Table 5 


Case ID 


Description 


39 


15001 Class 1 


40 


15000 Class 2 


41 


15000 Class 3 


42 


14997 Class 4 


43 


14999 Class 5 


44 


10410 Class 6 


45 


85 Class 7 


46 


10343 Class 8 


47 


82 Class 9 


48 


9509 Class 10 


49 


15 Class 11 


50 


15 Class 12 


51 


10 Class 13 


52 


12 Class 14 


53 


2 Class 15 



• t 



52 



Table 5 


Case ID 


DescriDtioii 


54 


1 Class 16 


55 


2 Class 17 


56 


1 Class 18 


57 


2 Class 19 



Fig. 21A - Fig. 21D illustrate traffic cases 20, 29, 39, 48 (again showing contours 
of Contours of q, B, and m), with maximum admissible and minimum inadmissible 
numbers of connections, average rate lOkbit/s. Fig. 21 A shows case 20 (constant rate, 
5 15000 connections); Fig. 21B shows case 29 (always ON or OFF, Pon=0.01, 9508 
connections); Fig. 21C shows case 39 (constant rate, 15001 connections); Fig. 21D 
shows case 48 (always ON or OFF, Pon=0,01, 9509 connections). 

Fig. 22A - Fig. 22D illustrate traffic cases 21, 22, 40, 41 (again showing contours 
of q, B, and m) for maximum admissible and minimum inadmissible numbers of 

10 cormections, average rate lOkbit/s, average period 10 ms. Fig. 22A shows case 21 ( 
ON-OFF periodic, Pon=0.01, 14999 connections); Fig. 22B shows case 22 (discrete 
periodic, 14999 connections); Fig. 22C shows case 40 (ON-OFF periodic, Pon=0 01> 
15000 coimections); Fig. 22D shows case 24 (Poisson, 15000 connections). In Fig. 22 
there is no difference between the ON-OFF periodic and discrete periodic models of 

15 approximately the same type of connections. This is because the burst size 100 bit is 
much smaller than the buffer size 3Mbit. 

Fig. 23A - Fig. 23D illustrate traffic cases 23, 24, 42, 43 (again showing contours 
of q, B, and m) for maximum admissible and minimum inadmissible numbers of 
connections, average rate lOkbit/s, average period 10 ms. Fig. 23A shows case 23 (ON- 
20 OFF Markov, Pon=0.01, 14996 connections); Fig. 23B shows case 24 (Poisson, 14998 
connections); Fig. 23C shows case 42 (ON-OFF Markov, Pon=0-01, 14997 
connections); Fig. 24D shows case 43 (Poisson, 14999 connections). 
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Fig. 24A - Fig. 24D illustrate traffic cases 25, 26, 44, 45 (again showing contoxirs 
of q, B, and m) for maximum admissible and minimimi inadmissible numbers of 
connections, average rate lOkbit/s, average period 100 s. Fig. 24A shows case 25 ( ON- 
OFF periodic, Pon=0.01, 10409 connections); Fig. 24B shows case 26 (discrete 
5 periodic, 84 connections); Fig. 24C shows case 44 (ON-OFF periodic, Pon=0.01, 10410 
connections); Fig. 24D shows case 45 (discrete period, 85 connections). In Fig. 24 
there is a great difference between the ON-OFF periodic and discrete periodic models 
of approximately the same type of connections. This is because the burst size 1Mbit is 
of the same order of magnitude as the buffer size 3 Mbit. 

10 

Fig. 25A - Fig. 25D illustrate traffic cases 27, 28, 46, 47 (again for contours of q, 
B, and m) for maximum admissible and minimum inadmissible numbers of connections, 
average rate lOkbit/s, average period 100 s. Fig. 25A shows case 27 (ON-OFF Markov, 
PoN=0.0l, 10342 connections); Fig. 25B shows case 28 (Poisson, 81 connections); Fig. 
15 25C shows case 46 (ON-OFF Markov, Pon=0.01, 10343 connections); Fig. 25D shows 
case 47 (Poisson, 82 connections). 

Fig. 26A - Fig. 26D illustrate traffic cases 30, 31, 49, 50 (again showing 
contours of q, B, and m) for maximum admissible and minimum inadmissible numbers 
of coimections, average rate lOkbit/s, average period 10 ms. Fig. 26A shows case 30 
20 (ON-OFF periodic, Pon=0.01, 14 connections); Fig. 26B shows case 31 (discrete 
periodic, 14 connections); Fig. 26C shows case 49 (ON-OFF periodic, Pon=0.01, 15 
connections); Fig. 26D shows case 50 (discrete period, 15 connections). 

Fig. 27 A - Fig. 27B illustrate traffic cases 32, 33, 51, 52 (again for contours of q, 
B, and m) for maximum admissible and minimum inadmissible numbers of connections, 
25 average rate lOMbit/s, average period 10 ms. Fig. 27A shows case 32 (ON-OFF 

Markov, Pon=0.01, 9 connections); Fig. 27B shows case 33 (Poisson, 11 connections); 
Fig. 27C shows case 51 (ON-OFF Markov, Pon=0.01, 10 connections); Fig. 27D shows 
case 52 (Poisson, 12 connections 
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Fig. 28A - Fig. 28D illustrate traffic cases 34, 35, 53, 54 (showing contours of q, 
B, and m) for maximum admissible and minimum inadmissible numbers of connections, 
average rate lOMbit/s, average period 100 s. Fig. 28A shows case 34 (ON-OFF 
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periodic, Pon=0 01, 1 connection); Fig. 28B shows case 35 (discrete periodic, 0 
connections); Fig. 28C shows case 53 (ON-OFF periodic, Pon=0.01, 2 connections); 
Fig. 28D shows case 54 (discrete period, 1 connection; same as case 16; no B <=3 Mbit 
and no m<=0 found) 

5 Fig. 29A - Fig. 29D illustrate traffic cases 36, 37, 55, 56 (showing contours of q, 

B, and m) for maximum admissible and minimum inadmissible numbers of connections, 
average rate lOMbit/s, average period 100 s. Fig. 29A shows case 36 (ON-OFF 
Markov, Pon=0.01, 1 connection; same as case 17); Fig. 29B shows case 37 (Poisson, 0 
connections); Fig. 29C shows case 55 (ON-OFF Markov, Pon=0.1, 2 connections); Fig. 

10 29D shows case 56 (Poisson; same as case 18). 



Fig. 30A - Fig. 30B illustrate traffic cases 38, 57 (showing contours of q, B, and 
m) for maximum admissible and ndnimum inadmissible numbers of connections, 
average rate lOMbit/s, always ON or OFF, Pon=0.1 Fig. 30A shows case 38, 1 
connection; same as case 19); Fig. 30B shows case 57 (2 connections). 

15 While the invention has been described in connection with what is presently 

considered to be the most practical and preferred embodiment, it is to be understood 
that the invention is not to be limited to the disclosed embodiment, but on the contrary, 
is intended to cover various modifications and equivalent arrangements included within 
the spirit and scope of the appended claims. 
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